%Tr - Receding contact angle
%Ta - Advancing contact angle
Tr_GO=[25.28 31.28 26.65 40.03]; %Tip 14 18 19 24
Ta_GO=[34.97 47.38 38.33 43.89];

Tr_1rGO=[55.16 36.11 45.10 43.67]; %Tip 14 18 19 24
Ta_1rGO=[59.12 40.89 54.48 47.10];

Tr_10rGO=[49.57 60.46 53.61]; %Tip 14 18 19
Ta_10rGO=[52.19 66.02 55.91];

x=[0.9 1.9 2.9];
Tr=[mean(Tr_GO) mean(Tr_1rGO) mean(Tr_10rGO)];
err_Tr=[std(Tr_GO) std(Tr_1rGO) std(Tr_10rGO)];
plot(x,Tr,'.r-','MarkerSize',15,'Color','r','LineWidth',0.7)
%errorbar(x,Tr,err_Tr,'linewidth',1)

hold on
% Tip 14
x=[0.85 1.85 2.85];
Tr=[25.28 55.16 49.57];
plot(x,Tr,'*r','MarkerSize',10)

hold on
% Tip 18
x=[0.85 1.85 2.85];
Tr=[31.28 36.11 60.46];
plot(x,Tr,'xr','MarkerSize',10)

hold on
% Tip 19
x=[0.95 1.95 2.95];
Tr=[26.65 45.10 53.61];
plot(x,Tr,'or','MarkerSize',8)

hold on
% Tip 24
x=[0.95 1.95];
Tr=[40.03 43.67];
plot(x,Tr,'sr','MarkerSize',10)

hold on
x=[1.1 2.1 3.1];
Ta=[mean(Ta_GO) mean(Ta_1rGO) mean(Ta_10rGO)];
err_Ta=[std(Ta_GO) std(Ta_1rGO) std(Ta_10rGO)];
plot(x,Ta,'.b-','MarkerSize',15,'Color','b','LineWidth',0.7)
%errorbar(x,Ta,err_Ta,'linewidth',1)
xticks([0 1 2 3 4])
xticklabels({' ','0','1','10',' '})
ylabel('Contact angle (Deg)')
xlabel('Reduction time (min)')
%legend('Receding','Advancing')

hold on
% Tip 14
x=[1.05 2.05 3.05];
Ta=[34.97 59.12 52.19];
plot(x,Ta,'*b','MarkerSize',10)

hold on
% Tip 18
x=[1.05 2.05 3.05];
Ta=[47.38 40.89 66.02];
plot(x,Ta,'xb','MarkerSize',10)

hold on
% Tip 19
x=[1.15 2.15 3.15];
Ta=[38.33 54.48 55.91];
plot(x,Ta,'ob','MarkerSize',8)

hold on
% Tip 24
x=[1.15 2.15];
Ta=[43.89 47.10];
plot(x,Ta,'sb','MarkerSize',10)